The phase diagram of neutral quark matter 
with pseudoscalar condensates in the color-flavor locked phase 
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We calculate the phase diagram of color and electrically neutral quark matter in weak equilibrium 
using the NJL model at high baryon densities. We extend previous analyses by taking into account 
the possibility of formation of pseudoscalar meson condensates in the color-flavor locked (CFL) 
phase. We find that the kaon condensed phase (CFL-K ) is preferred over the CFL phase at low 
temperatures, ft is found that the CFL-K phase contains no gapless modes. The results we have 
obtained should be of relevance for the modeling of compact stars. 



I. INTRODUCTION 

At sufficiently high baryon densities, massless three- 
flavor QCD is in a color-flavor locked (CFL) phase Q . In 
the CFL phase the original SU(3) L x SU(3)ij x SU(3) C x 
U(1)b symmetry of QCD is broken to SU(3)y +c due to 
the formation of scalar diquark condensates. As a conse- 
quence all eight gluons become massive by the Anderson- 
Higgs mechanism. Furthermore the breakdown of chiral 
symmetry leads to the emergence of eight massless pseu- 
doscalar excitations, similar to what happens in QCD at 
low baryon densities. Casalbuoni and Gatto have con- 
structed a chiral effective theory for the CFL phase Q 
that describes these massless excitations. In such a chiral 
effective theory the pseudoscalar excitations 7r a are con- 
tained in a flavor octet chiral field S = exp(i7r a A a // 7r ). 
Here / w is the pion decay constant and the Cell-Mann 
matrices A a are normalized as tr A a Af, = 25 a t„ where a = 
1 . . .8. The field £ transforms under SU(3) L x SU(3) fl 
as £ -► Ul^U^. 

In the CFL phase with massless quarks the vacuum 
expectation value of the S field is equal to 1. How- 
ever, finite quark masses and chemical potentials give 
rise to a potential for E [1, 0, 0- Depending on the 
values of the chemical potentials and the current quark 
masses m u , rrid and m s the new ground state can have 
S ^ 1, which arises via an axial flavor transformation on 
the CFL ground state. In such a case one of the pseu- 
doscalar meson fields obtains a vacuum expectation value 
P| . Depending on the fields involved one speaks of pion 
or kaon condensation giving rise to the CFL-tt*, CFL-K 
and CFL-K± phases 

Using the chiral effective theory it was found that 
pseudoscalar condensation can occur if the chemical po- 
tential of the mesonic excitation exceeds the mass of the 
corresponding excitation. For the case m u = rrid the 
masses of the pionic and kaonic excitations (not to be 
confused with the masses and mx of the pion and 
kaon mesons) in the chiral effective theory for the CFL 



phase are respectively found to be Q M 2 = 2am u m s and 
= a(m u + ms)m u . At asymptotically high densities 
a — 3A 2 /(7r 2 / 2 ) (3, where A is the CFL gap and the pion 
decay constant is given by / 2 = (21 — 8 log 2)^ 2 /(367r 2 ). 
Here (j, denotes one third of the baryon chemical poten- 
tial. The chemical potentials of the pseudoscalar mesons 
are @ /v+ = HQ, 



Hk+ — HQ 
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where fj,Q is the chemical potential for electric charge. 
The chemical potentials of the antiparticles are the op- 
posite of those above. The 7T° does not condense since 
its chemical potential vanishes [||. The rj and the rf ex- 
citation can however obtain a vacuum expectation value 
but will not be taken into account in this article. 
The parameters of the chiral effective theory ana- 
lysis have to be matched to QCD, which only can be done 
accurately at extremely high baryon densities where per- 
turbative QCD is applicable. However, it is interesting 
to know how matter behaves at non-asymptotic baryon 
densities, such as inside a compact star. Therefore it 
would be worthwhile to have an alternative analysis of 
pseudoscalar condensation in the color superconducting 
phases which could be reliable at intermediate baryon 
densities. The model we use for that purpose here is the 
Nambu-Jona-Lasinio (NJL) model (see y| for a review). 
The NJL model has been applied to study the QCD phase 
diagram as a function of baryon chemical potential and 
temperature under compact star constraints 0, lid ITl) 
(see e.g. 0] for a review of compact stars). In this ar- 
ticle we will extend the NJL model analysis of the phase 
diagram of Refs. by including the possibility of 

pseudoscalar condensation in the color-superconducting 
phases. In this way we are able to obtain a more detailed 
phase diagram for quark matter under compact star con- 
straints, that is electrically and color neutral matter in 
weak equilibrium. Earlier investigations of pseudoscalar 
condensation in the CFL phase u sing the NJL model 
have been performed in Refs. [H [l||. The phase di- 
agram of quark matter under compact star constraints 
was however not calculated in Refs. [Tol IT^ |. Reason- 
able agreement between the chiral effective theory and 



the NJL model approach to pseudoscalar condensation 
in the CFL phase was found in [l3l IhH ] . This we will use 
as a justification to apply some of the outcomes of the 
chiral effective theory analysis. 

This article is organized as follows. In Sec. II we dis- 
cuss shortly the NJL model and the parameter choices. 
The details of the calculation are given in Sec. III. The 
reader only interested in the main results may go directly 
to Sec. IV in which the phase diagram of quark matter 
under compact star constraints is presented. We summa- 
rize and conclude in Sec. V. 



II. THE NJL MODEL 

In the NJL model, one treats the interaction between 
the quarks as a point-like quark color current-current in- 
teraction. We will use the following Lagrangian density 

C = i> {i^dft - m + p,j ) ip + Cq q + C qq + £h, (2) 

where the quark-antiquark term, C qq , the diquark inter- 
action term, C qq , and the 't Hooft interaction term £h 
are defined below. We have suppressed the color (c), fla- 
vor (f), and Dirac (d) indices of the fermion fields ip for 
notational simplicity. The mass matrix m is diagonal and 
contains the current quark masses. The matrix fi is also 
diagonal in flavor space and contains the quark chemi- 
cal potentials. In order to allow for electric and color 
neutralization 0, 0, 0, we write 

p, = nl c <8 lf+fiQl c <8 Q + fah + ^8*8 01/, (3) 

where Q = diag(2/3, — 1/3, — 1/3). When we later on in 
Sec. Ill we introduce a electron (e) and muon (fj,) gas, this 
choice will automatically enforce the weak equilibrium 
conditions, (id,s = M« + Me,f»- The color matrices corre- 
sponding to the color chemical potentials /13 and are 
t 3 = diag(l,-l,0) and t s = diag(l, 1, -2)/ v / 3. We use 

the metric g^ v — diag(H ) and the standard Dirac 

representation for the 7-matrices. The quark-antiquark 
interaction part of the Lagrangian density is 



Co 



G 



(lpX a 1p) + (lp\ a i~f 5 lp) 2 



(4) 



Here the matrices A a denote the 9 generators of U(3) and 
act in flavor space. They are normalized as TrA Q Ab = 
25 a b- To remind the reader, the antisymmetric flavor 
matrices A2, A5 and A7 couple up to down, up to strange 
and down to strange quarks, respectively. The diquark 
interaction term of the Lagrangian density is given by 



C qq = -G (^t A \ B Ci l5 ^ T ) (iP T t A \ B Ci l5 if) 
+-G (?pt A \ B Ci> T ) (iP T t A X B Cip) , 



(5) 



where A, B S {2, 5, 7} since only the interaction in the 
color and flavor antisymmetric triplet channel is attrac- 
tive. The matrices t a are the generators of U(3) and act 



in color space. Their normalization is tri a ifc = 2S a b- The 
charge conjugate of a field ip is denoted by tp c — Cip T 
where C = 17270- The coupling strength 3G/4 of the di- 
quark interaction is fixed by the Fierz transformation (see 
for example Ref. |8(). However, some authors discuss the 
NJL model with a different diquark coupling constant (in 
Ref. @ it was found that changing this coupling strength 
to G can have a non-negligible effect on the phase dia- 
gram) . Finally, the 't Hooft interaction which models an 
additional breakdown of the axial U(l) symmetry due to 
instantons is given by 0] 

£ H = -#[det / $(l-7 5 )^ + det / ^(l+7 5 )^] . (6) 

The results that will be presented below are obtained 
with the following choice of parameters 0, 0] 

m u = m d = 5.5 MeV, m s = 140.7 MeV 
G= 1.835/A 2 , K = 12.36/A 5 , A = 602.3 MeV. (7) 

This choice of paramete rs g ives the following parame- 
ters the realistic values |19| = 135.0 MeV, mx = 
497.7 MeV, m v = 514.8 MeV, ny = 957.8 MeV, and 
fir — 92.4 MeV. Moreover, it allows us to compare our 
results to the phase diagram obtained in 0, where only 
the scalar diquark condensates were taken into account. 



III. CALCULATION OF THE PHASE 
DIAGRAM 

To obtain the phase diagram of the NJL model, we 
introduce the chiral condensates 



<j u = (uu) , a d = (dd) , 



cr, = (ss) 



(8) 



and the scalar and pseudoscalar diquark condensates 



t A X B C l5 ij), A p AB = —G(ip T t A X B Cip). 

(9) 

We do not take into account the possibility of pion con- 
densation or charged kaon condensation outside the CFL 
phase (e.g. the appearance of (uvy^d) and (uij^s), con- 
densates). Using the results of |2(| this is in principle 
possible if \/iq\ > or \hq\ > tor- Phase diagrams 
with these condensates have been investigated in Refs. 
Ullliil. From Ref. [|| we conclude that elec- 
tric neutralization gives values of \hq\ which are smaller 
than mx making the occurrence of this type of kaon con- 
densate unlikely. Electric neutralization allows for values 
\/j,q\ which are larger than m w , but from the analysis of 
Ref. [2^| it follows that for matter in weak equilibrium the 
phase with the pion condensate is not competing against 
color superconducting phases. This makes it in our opin- 
ion improbable that such a condensate arises in neutral 
quark matter at high baryon densities. 

In the CFL phase with massless quarks the nonzero 



condensates are A, 



In the more realistic 



situation of unequal quark masses and unequal chemical 



2 



potentials the same condensates appear, but then they 
are no longer equal to each other. As was explained in 
the introduction the state with pseudoscalar condensa- 
tion in the CFL phase can be obtained by an axial fla- 
vor transformation on the CFL state. In such a case 
ip — ► exp(i75# a A Q /2)?/> and the scalar diquark conden- 
sates are partially rotated into pseudoscalar diquark con- 
densates [laH. For a K° mode the diquark condensates 
transform as 



i o5 



A 



77 



cos(0/2)A* 2 , 
cos(0/2)A§ 5 , 
•A? 



A£ 5 
Af 2 



x 77i 



ism(0/2)(6 6 ~i9 7 )A s 22 , 
ism(6/2)(§e + ie 7 )A s 55 , 

(10) 



where 8 — \JQ\ + 7 and 8 a = 6 a /9. In the same way A| 2 
and A 77 are partially rotated into A 27 and A 72 under a 
K* transformation, while ir^ transformations rotate A| 5 
and A 77 into Ag 7 and A 75 [J^ . According to the analysis 
in the chiral effective theory and confirmed numerically 
by Buballa in the NJL model, the appearance of 
a pseudoscalar diquark condensate lowers the effective 
potential by Q 



/^ 2 (l-cos#) 2 /2, 



(11) 



where i = K°, K* or tt^ and cos 9 = Mf /fj, 2 . 

So the CFL-K phase is characterized by the CFL 
condensates plus two non- vanishing pseudoscalar diquark 
condensates A 25 and A^ 2 . Similarly, CFL-K* and CFL- 
7r ± phases exist. Next to the CFL phase, different 
color superconducting phases are known, like the 2SC 
|2g (two-flavor color superconducting) phase (A| 2 ^ 
0,A§ 5 = Af 7 = 0) and the uSC phase (A| 2 ^ 
0, Ag 5 ^ 0, A 77 = 0). In the phase diagram we will 
present, we also encounter a phase with a A 22 conden- 
sate together with a non- vanishing Af 2 . This phase we 
will call p2SC. A nonzero value of Ag 2 can be obtained by 
an axial color transformation on the 2SC state. For con- 
venience we have summarized the possible combinations 
in Table IJ 

To obtain the phase diagram of the NJL model 
under compact star constraints with the possibility of 
pseudoscalar condensation in the color-superconducting 
phases we will proceed as follows. In order to allow 
for electric neutralization, we introduce a free electron 
and muon gas and calculate the finite temperature (T) 
mean field effective potential which yields (see for exam- 
pie iH) 



V = 2G(a 2 u +a 2 +a 2 ) 



IA 



AB\ 



+ \Kb\ 



3G 



T . f d 3 p 
- 4K<T u a d cr s - — ^2 / 7 — 3 logdetS*" 1 



(27T)- 



(12) 



where Q n = (2n + X)irT and E t = y/p 2 + mf. Weak 
equilibrium is automatically satisfied in this way (for 
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TABLE I: Definition of different phases. A or a x denotes 
respectively a vanishing or non- vanishing condensate. If at 
least one of the quasi-particles has a gapless mode one writes 
a 'g' in front of the name of the phase, e.g. g2SC, gCFL, etc. 
If fiQ > 0, one speaks of CFL-K + and CFL-7r + , if /iq < 
of CFL-K" and CFL-rr". The five phases above the horizon- 
tal line appear in the phase diagram of neutral quark mat- 
ter, Fig. |3 The other phases are possibilities (among others) 
which we have taken into account, but did not arise in the 
phase diagram. 



example one can read off that fi s = fi u + /x e ). Here 
we have assumed that neutrinos leave the star immedi- 
ately. An analysis with neutrino trapping, as done for 
the phase diagram with scalar diquark condensates in 
Ref. [25| . is left for future work. The electron and muon 
mass have the following values, m e = 0.511 MeV and 
to m = 105.66 MeV. The 72 x 72 matrix S -1 contains the 
inverse Nambu-Gorkov fermion propagator in the mean 
field approximation 



l c <g> V + ft <g> 70 t A ®\ B 
t A ®\ B ®$ 2 AB l c ®£>-A®7o 



(13) 



Here 



and M = diag(mo M — 4Ger„ + 2Ka d a, 
2Ka s cr u , m 0s - AGa s + 2Ka u a d ) 
and & AB are equal to 



(14) 

m od - AGa d + 
The matrices & AB 



$1 



AB 
2 

AB 



= A^ B 7 5 4 
= -A^75 



A^ B U 



(15) 
(16) 



The matrix 1 is the identity matrix in color (c), flavor 
(/), or Dirac (d) space. 

To calculate the effective potential in an efficient 
way, one can multiply the matrix S" -1 with diag(l c £g) 
l/<g>7o,l c ( S>l/ ( S>7o) which leaves the determinant in- 
variant. In this way, one obtains a new matrix R with 
ipo's on the diagonal. We can write R — ipol + A, where 
A is a Hermitian matrix independent of po- By deter- 
mining the eigenvalues of the matrix A, one can recon- 
struct the determinant of S^ 1 for all values of po which 



3 



is nl=i(^» + Wo)- After summing over Matsubara fre- 
quencies, one finds 



T^logdetS- 1 = 

PO=^n 



72 

E 

i=l 



^ + Tlo g 



l + e-l^l/T 



(17) 

All that remains in order to determine the effective 
potential is to perform the integration over three- 
momentum p up to the ultraviolet momentum cutoff A 
numerically. 

In order to enforce electric and color charge neutral- 
ity the following equations have to be satisfied 



dV 



= 0, n 3 ,8 



dV 



= 



3.8 



(18) 



where uq is the electric charge density, and 77.3 = n r — n g , 
and ns = (n r + n g — 2n; ) )/v / 2. 

The values of the condensates and the phase diagram 
are determined by minimizing the effective potential V 
with respect to the set of possible condensates. Together 
with the electric and charge neutrality constraints this 
implies that we have to solve the following equation 



dV 

dx. 



= 0, 



(19) 



where x = {a u , a d , a s , A S AB A P AB , fi Q , ^ 8 }. 

We will allow for the chiral condensates and the pos- 
sible diquark condensates denoted in Table [I] The analy- 
sis in the effective theory shows that the different phases 
with pseudoscalar condensation are separated by a first- 
order transition 6J. Therefore it is reasonable to mini- 
mize the effective potential with respect to the different 
possible sets of pseudoscalar diquark condensates sepa- 
rately. Thereafter one can compare the minimal values 
of the effective potential in order to determine the phase. 
This speeds up the numerical calculation. 

Using the diagonal color transformations and the 
U(l) flavor transformations it is always possible to choose 
the scalar diquark condensates to be real, the phases 
of the pseudoscalar condensates however cannot be ro- 
tated away. We find that the effective potential of the 
NJL model with real CFL condensates and two pseu- 
doscalar diquark condensates A?^ and Af? 2 is invari- 
ant under the transformation — > e'^A^ together 
with Ag 2 — > e~^A5 2 (for the other possible sets of 
diquark condensates, the same transformation rule of 
course holds). This shows that one can always rotate 
one diquark condensate to be real. From Eq. I|10[l it fol- 
lows that in the minimum of the effective potential the 
other diquark condensate is then real as well. We have 
also checked this numerically, by allowing for a complex 
pseudoscalar diquark condensate. Hence we can take the 
three scalar and the two pseudoscalar diquark conden- 
sates to be real. This speeds up the calculation as well. 

To evaluate the derivatives of the effective potential 



(see also [TT| ) we have used the following formula 



d 



T dx~ E logdetS-V) 



72 , 



e |A,|/T + i 



sgn(A), (20) 



where hi = (U> dA/dxjU)u. Here Ai are again the eigen- 
values of A, and U is a unitary matrix which contains 
in the i-th column the normalized eigenvector of A with 
eigenvalue Xi . To obtain the complete derivative of the 
effective potential one still has to integrate Eq. I|20|) over 
momenta up to the cutoff A. The advantage of using 
Eq. H20JI to evaluate the derivatives is threefold: firstly it 
is more accurate than the finite difference method, sec- 
ondly it is also much faster since one had to diagonalizc 
the matrix A anyway in order to obtain the value of the 
effective potential, and finally additional derivatives can 
be computed without much extra numerical work. 

The speed of the calculation of the effective potential 
and its derivatives depends heavily on how fast one can 
compute the eigenvalues. There are several ways to speed 
up the calculation. Firstly, the determinant of A does not 
depend on the direction of p. Therefore, one can choose p 
to lie in the z-direction. Together with the choice of the 
non- vanishing condensates mentioned above, this implies 
that A becomes a real symmetric matrix, which simpli- 
fies the calculation of the eigenvalues. Secondly, one can 
interchange rows and columns of A without changing its 
determinant. By doing so, one can bring A in a block- 
diagonal form. One can then determine the eigenval- 
ues of the blocks separately which is significantly faster 
since the time needed to compute eigenvalues numeri- 
cally scales cubically with the dimension of the matrix. 
In the most general case with two pseudoscalar diquark 
condensates, one can always reduce the problem to one 
20 x 20 matrix and three 8x8 matrices. 

We determined the eigenvalues using LAPACK 
routines |28j| . The numerical integration over three- 
momentum p up to the cutoff was done using a Gauss- 
Legendre quadrature with 2 x 32 points at T > 15 MeV. 
For lower temperatures the Bose-Einstein distribution 
functions start to behave like step-functions, therefore 
the Gauss-Legendre quadrature, which is based on poly- 
nomial interpolation, becomes less accurate. For that 
reason we used Simpson integration with 256 points at 
lower temperatures. We have checked the accuracy of the 
integration procedures by doubling the number of inte- 
gration points by 2. The condensates were determined 
by solving Eqs. I|19fl an d i|18|) using a multidimensional 
rootfinding algorithm of the GSL package 29] with differ- 
ent initial conditions. The ground state was determined 
by picking the solution that has the lowest value of the 
effective potential. After the ground state was found we 
investigated the zeros of |Aj| as a function of momentum. 
If I Ai| has a zero, one of the quasi-particles has a gap- 
less mode. In that case one speaks of the gCFL, g2SC, 
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or guSC phase |30(. The problem of gapless phases are 
its chromomagnetic instabilities . Therefore such a 
phase cannot be the phase with the lowest free energy. 



IV. RESULTS 

In this section we will discuss the phase diagram of 
quark matter under compact star constraints, that is elec- 
trically and color neutral matter in weak equilibrium. We 
will discuss two possible checks of our results. The first 
is to compare with the results of the chiral effective the- 
ory. The other is to compare with the phase diagram of 
Ref. in which no pseudoscalar color-superconducting 
phases were taken into account. 

We first present the phase diagram in the /x_b-/xq 
plane (so no electric neutralization) at T = 0. This al- 
lows for comparison with the chiral effective theory cal- 
culations of |6|. Thereafter we will discuss the phase 
diagram of neutral quark matter. 



A. The hb-1-iq phase diagram at T = 

To predict the shape of the ^b~^q phase diagram 
one can use Eqs. (JTJ and The resulting phase dia- 

gram is somewhat different from the analysis in the chi- 
ral effective theory of Ref. [g , because here electromag- 
netic corrections to the masses of the charged mesonic 
excitations have not been taken into account. Firstly, 
if /iq is positive, /kk+ > Mk° and /ik+ > Mtt+- Since 
the energy gain of meson condensation is proportional to 
nl-2Mf+Mf/nf and m > Mi, it follows that the CFL- 
K + phase is favored over the CFL-K°/7r + phase for posi- 
tive values of /iq. For negative values of /iq the CFL-K 
phase is favored over the CFL-K" phase, because then 
/!k- < Hk°- So there has to be a phase transition from 
the CFL-K phase to the CFL-K+ phase at = 0. Sec- 
ondly, if fiQ is negative enough, pion condensation will be 
favored over neutral kaon condensation. Since Mr, 
and /xko are all proportional to 1/fj, the transition line 
should be proportional to as well at asymptotically 
high densities. 

These conclusions are in agreement with the calcu- 
lation of the phase diagram presented in Fig. ^ For 
/iq > we find the CFL-K+ phase. There is a first-order 
phase transition at fiq = 0. At negative /iq we encounter 
first the CFL-K phase and then the CFL-7T - phase. If 
the baryon chemical potential is increased, the transition 
from the CFL-K to the CFL-7r~ phase occurs at smaller 
values of /iq, in agreement with the predictions. 

We want to stress that, in order to calculate this 
diagram, we have used that chiral effective theory analy- 
sis [6j shows that the phase boundaries between the the 
phases with pseudoscalar condensation are first-order. 
Therefore, in order to simplify the calculation, we have 
assumed that this holds in the NJL model as well. 
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FIG. 1: The phase diagram of color neutral quark matter in 
weak equilibrium, as a function of baryon chemical potential 
and the electric charge chemical potential for T = 0. The 
solid lines denote first-order phase transitions. 



The phase diagram of color neutral matter in weak 
equilibrium as presented in Fig. ^ does not necessarily 
contain all phases which are present in the phase diagram 
of matter that is also electrically neutral. This is because 
the electric neutrality constraint can force the electrically 
neutral system to be in a state which is a meta-stable 
state of color neutral matter in weak equilibrium. 



B. The phase diagram of neutral quark matter 

In Fig. El we have displayed our main result, the 
phase diagram of quark matter under compact star con- 
straints, that is color and electrically neutral matter in 
weak equilibrium. This phase diagram was obtained by 
minimizing the effective potential on a grid of 440 (in the 
fiB direction) by 300 (in the T direction) points. In this 
way we can determine the phase boundaries up to an ac- 
curacy of ±0.25 MeV in the direction and ±0.1 MeV 
in the T direction. After we minimized the effective po- 
tential on this grid we have checked the continuity of the 
minimum value of the effective potential. 

In Fig. |3| we have displayed the values of the con- 
densates and the chemical potentials in the minimum of 
the effective potential for T = MeV, T = 20 MeV and 
T = 40 MeV. Outside the CFL phase, excellent agree- 
ment is found with Q (their ^ and /Js differ by a factor 2 
from ours due to another normalization). Since we took 
the pseudoscalar condensates to be real, a characteristic 
feature can be read off from Eq. I|I0|I , namely that they 
should have opposite sign. This we have indeed found 
and can be seen in our results in Fig|3J 



The most important conclusion to be drawn from the 
phase diagram is that (for our realistic parameter choices) 
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FIG. 2: The phase diagram of electrically and color neutral quark matter in weak equilibrium, as a function of baryon chemical 
potential and temperature, first-order phase transitions are denoted by a solid line, second-order phase transitions by a dashed 
line. The dotted lines indicate the (dis)appearance of gapless modes in the quasi particle excitation spectrum. Critical endpoints 
are denoted by a black dot. The xSB and the NQ labels denote the chirally broken and neutral chirally symmetric quark matter 
phases respectively. The definition of the color- superconducting phases can be found in Table [I] 



electrically and color neutral quark matter in weak equi- 
librium is at high baryon densities and low temperatures 
in the CFL-K phase. We find that this phase contains 
no gapless modes, which implies that it does not contain 
chromomagnetic instabilities. 

The CFL-K is not always favored over the CFL 
phase. If at high baryon densities the temperature is 
raised the CFL-K phase goes via a second-order phase 
transition to the CFL phase. At lower baryon densities 
we find a critical point where the transition to the CFL 
phase becomes first-order. 

In electrically neutral quark matter, there should al- 
ways be an excess of electrons and muons to compensate 
for the lower density of strange quarks compared to up 
and down quarks. Therefore /iq has to be necessarily 
negative (this can be seen in Fig. [31 as well), which im- 
plies that only the CFL-K and the CFL-7T phase could 
in principle be found in electrically neutral quark matter. 
We find the CFL-K phase, but we do not encounter the 
CFL-7r~ phase in the phase diagram. 

Another striking feature of the phase diagram is 
the occurrence of the p2SC phase, a two-flavor color su- 
per conductor with a non- vanishing pseudoscalar diquark 
condensate, A| 2 . We find that the difference in free en- 
ergy between the p2SC phase and the 2SC phase is excep- 
tionally small. While usually a difference in the minimal 
value of the effective potential between two competing 
phases is in the order of 0.1 — 1 MeV/fm 3 , in this case 
it is only of order 0.1 keV/fm 3 , see Fig. 0] Therefore it 
could well be that other parameter choices will let the 
p2SC phase disappear, but nevertheless it is very inter- 
esting that its free energy is so close to the ordinary 2SC 



phase. 

We have displayed the values of the condensates 
of the p2SC phase in Fig. again for T = 30 MeV. 
From this figure one can see that upon increasing the 
baryon chemical potential, one enters the p2SC phase via 
a second-order phase transition. The dotted line in the 
figure indicates the A| 2 condensate of the meta-stable 
2SC phase. One can see from Fig. [S] that the A| 2 con- 
densate has a relatively small value compared to the A| 2 
condensate. We find that the square root of the sum of 
the values of the A| 2 and the A 25 condensates is equal 
to value of the meta-stable A| 2 condensate. This shows 
that the p2SC phase can be obtained from the 2SC phase 
by an axial color transformation. 

In our definition the p2SC phase contains a non- 
vanishing A| 2 and a Ag 2 condensate. Using the off- 
diagonal color rotations, one can rotate a Ag 2 in a A^ 2 
condensate. We have checked that taking into account 
a A| 2 and A^ 2 condensate leads to the same phase dia- 
gram. We have also investigated whether other possible 
combination of A 22 , A?? 2 and A^ 2 lead to a lower free 
energy, but we did not find such solutions with lower free 
energy. 

We have also sought other possible combinations of 
two pseudoscalar diquark condensates and three scalar 
diquark condensates (of which one or more could poten- 
tially vanish), but we did not find any of them. 

Comparing our results to the phase diagram of 
(where the CFL-K phase was not considered) we have 
good agreement outside the CFL-K and the p2SC phase. 
The phase boundaries of the 2SC phase coincide. Also 
we find the boundaries of the g2SC and the guSC gapless 



6 



T = MeV 




300 
200 
100 

10 
5 

-5 
-10 
-15 
-20 



350 400 450 500 550 



M3 



MS 

—I I— 



350 400 450 500 550 




350 400 450 500 550 



60 
40 
20 

-20 
-40 
-60 
-80 
-100 



1 1 

. — MQ 


1 1 


-- A£ 5 










/ 











350 400 450 500 550 

Ms/3 (MeV) 



600 
500 
400 
300 
200 
100 

10 



T = 20 MeV 



— 1 1 


— 1 1 

M u — 




M d - ■ 




... M s 















600 
500 
400 
300 
200 
100 




T = 40 MeV 



— 1 1 


— 1 1 

M u — 




M d - ■ 




M s 















350 400 450 500 550 



350 400 450 500 550 



-10 









1 






/ 

/ 
1 
1 
1 




0.5 



v. 
/ •• 








\ A ; 




\ 






-- M3 

A*8 

i i 


i 


_i 


-0.5 
i 


M3 

Ms 



350 400 450 500 550 



350 400 450 500 550 





350 400 450 500 550 




350 400 450 500 550 
Ms/3 (MeV) 



350 400 450 500 550 

Ms/3 (MeV) 



FIG. 3: Constituent quark masses, chemical potentials and condensates for electrically neutral quark matter in weak equilibrium. 
The constituent quark masses, chemical potentials and condensates are displayed as a function of baryon chemical potential, 
for T = MeV (left-hand side), T — 20 MeV (middle panels) and T = 40 MeV (right-hand side). 
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phases at the same place. At higher temperatures this 
also holds for the CFL and the gCFL phase boundaries 
(at lower temperatures we find the CFL-K phase). In- 
spection shows that at zero temperature at fj, = 443 MeV 
there is a first-order transition from the neutral quark 
matter phase to the CFL-K phase. This transition 
oint is in excellent agreement with the calculations of 
where at the same point a transition to the CFL 
phase was found. A difference is that we find that 
there is a very small window for the uSC phase, even 
for T < 5 MeV. Furthermore, our critical endpoints lie 
at somewhat higher temperatures. 

It must be said that other sensible choices of parame- 
ters can lead to different phase diagrams (see e.g. [9j for a 
comparison) . We however believe that the CFL-K phase 
should be relatively stable against changing parameters. 
This does probably not hold for the p2SC phase, because 
of the very small difference in free energy between the 
p2SC and the 2SC phase. One must also be aware that in 
reality also crystalline phases [33| and phases with color- 
spin locking (see e.g. |23) are possible. Also effects of 
non-local interactions could potentially modify the phase 



structure [34(. At low baryon densities, nuclear effects 
become important, which should also be taken into ac- 
count (see e.g. [3f| for a unified description of nuclear 
and quark matter in the NJL model). 

Let us finally indicate some astrophysical conse- 
quences of the phase diagram displayed in Fig. [2] com- 
pared to the phase diagram obtained in [t| in which no 
CFL-K phase was taken into account. In [jj the CFL- 
K° phase is replaced by a CFL and a gCFL phase. The 
difference in free energy between the CFL-K and the 
CFL phase is not large, so the new phase may not have 
a huge influence on the equation of state and hence the 
mass-radius relationships of a compact star. However, 
the neutrino cooling of a star is very sensible to the low- 
energy excitations [36j. In phases without gapless modes, 
like the CFL-K phase we found, the neutrino emissivity 
is exponentially suppressed compared to in phases with 
gapless modes [3(|. Hence replacing the gCFL by the 
CFL-K phase will reduce the cooling speed of a star 
with a quark core. 



V. CONCLUSIONS 

In this paper, we have calculated the \ib vs. T phase 
diagram of electrically and color neutral quark matter 
in weak equilibrium, at high baryon densities. It could 
be useful for understanding compact stars with a quark 
matter core. 

Our analysis has shown that neutral quark matter 
is in the CFL-K phase at high baryon densities and low 
temperatures, rather thane the CFL phase. Other CFL 
phases with pseudoscalar condensation were not found. 
This CFL-K phase is different from the ordinary CFL 
phase, since in the CFL-K phase two pseudoscalar di- 
quark condensates appear. The CFL-K phase we ob- 
tained, contains no gapless modes. This implies that 
there are no chromomagnetic instabilities in the CFL-K 
phase and that the neutrino emissivity is exponentially 
suppressed in the CFL-K phase. 

The CFL-K phase does not completely replace the 
CFL phase. At higher temperatures we found a first and 
a second-order phase transition to the CFL phase. 

Next to the CFL-K phase, we found at intermediate 
densities and temperatures a 2SC phase which contains 
a pseudoscalar diquark condensate, which we have called 
p2SC. This phase can be obtained by an axial color rota- 
tion on the 2SC phase. The difference in the free energy 
between the 2SC and the p2SC phase is very small. This 
implies that corrections may let this phase disappear. In 
any case it would be interesting to have a more detailed 
investigation of this phase. 

We have obtained good agreement with earlier cal- 
culations of the phase diagram were pseudoscalar con- 
densation was not taken into account and with results 
from the chiral effective theory. 
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